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Abstract 



The magnetic polarization induced by nonmagnetic impurities such as Zn 
in high Tc cuprate compounds is studied by the variational Monte Carlo sim- 
ulation. The variational wave function is constructed from the eigenstates 
obtained from Bogoliubov de Gennes mean field Hamiltonian for the two- 
dimensional t — J model. A Jastrow factor is introduced to account for the 
induced magnetic moment and the repulsion between holes and the impurity. 
A substantial energy gain is obtained by forming an antiferromagnetic polar- 
ization covering 4 or 5 lattice sites around the impurity. We also found the 
doping dependence for the induced magnetic moment consistent with experi- 
ments. 
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Recently a number of experiments, the neutron scattering^, nuclear-magnetic-resonance 
(NMR)i'i and scanning tunneling microscopy (STM)H, have been carried out to study 
the impurity effect on the electronic transport and magnetic properties in high T c cuprate 
compounds. These studies provide a detail information about the relationship between mag- 
netism and superconductivity in high T c cuprates. The nonmagnetic impurity Zn was found 
to suppress T c more strongly than magnetic impurity Ni, even though both replace Cu in the 
Cu0 2 pland. The amazingly accurate measurement of the local density of states (LDOS) 
by STMH also provides very different spectra for Zn and Ni. The spin dynamics studied by 
the neutron scattering experiments reveals that the low-energy spin fluctuations are strongly 
enhanced near the impurity and the magnetic excitation at the antiferromagnetic wave vec- 
tor (tt, tt) disappears with Zn doping in the underdoped regioni'i. It is interesting to find 
from the NMR and SQUID experiments that both the nonmagnetic Zn and the magnetic Ni 
impurities induce a local magnetic moment on Cu sites surrounding the impurity in the nor- 
mal state. The broadening of 63 Cu and 17 NMR lines has been attributed to a distribution 
of magnetic moments or a spatially inhomogeneous spin polarization extending over several 
lattice sites around the impurity. On the other hand some experiments^ found no evidences 
of the existence of local magnetic moments, at least in the optimum and overdoped samples. 
More careful theoretical and experimental efforts to exmine the magnetic polarization are 
needed to clarify this issue. 

So far most of the theoretical work has been based upon phenomenological BCS type 
models with emphasis on understanding of the LDOS. The observed nearly-zero-energy- 
resonance peak near Zn impurity was explained very early by Balatsky, Salkola and co- 
workers0Hll by assuming Zn to be an unitary impurity. Studie&§ based upon t — J type 
models have also successfully explained the LDOS. There are only few studied about the 
structure of magnetic polarization induced by the magnetic moment binded to the nonmag- 
netic impurity and the screening of this moment by other electrons. However in a recent 
paper0 Tsuchiura et al. use Gutzwiller approximation and the Bogoliubov-de Gennes (BdG) 
approach for the t — J model and they find no evidence of the existence of the local moments 



around the Zn impurity. They also concluded that the electron avoids the impurity instead 
of being binded to it. A much more careful examination of the effect of a non-mgnetic 
impurity in the t — J model is needed to resolve the controversy. 

Comparing with other phenomenological models, the t — J model has much stronger 
magnetic correlation and it may lead to a different picture about the magnetic polarization 
around the impurity. However, previous studies of the t-J model use the BdG approach with 
or without the Gutzwiller approximation and the no-doubly-occupied constraint imposed by 
the t — J model is only taken into account on the average or approximately. It very likely 
underestimates the antiferromagnetic correlation inherent in the t — J model. Another issue 
has not been addressed adequately before is the doping dependence of the induced magnetic 
moment. Very different results reported by NMR experimentsitE may be related to the 
doping dependence. 

In this paper we will impose the constraint rigorously by using the variational Monte 
Carlo approachil to study the effect of nonmagnetic Zn impurity on the ground state of 
the t — J model. The ground state trial wave function is first constructed by assuming 
d-RVB order parameters in the BdG approach. Then the variational wave function is shown 
to be greatly improved by adding a Jastrow factor to account for the strong magnetic 
correlation. We found a large energy gain by having an antiferromagnetic polarization 
around the impurity with size about 4 to 5 lattice sites as observed in 63 Cu NMR data! 
in the underdoped region. The significant suppression of the magnitude of the induced 
moment and its polarization size as doping increases to optimum doping is also consistent 
with experimental observations!!!. In addition, our result also provides a reason to explain 
the similarity between results! measured for Li + and Zn 2+ . Contrary to the work reported 
in Ref.(21) we show that electrons are always attracted to the impurity. But the effect gets 
weaker when number of holes increases. 

The model we consider is the dilute impurity limit of the two-dimensional t — J model. 
The interaction between impurities is neglected. Zn 2+ [3d 10 ] has total spin S = and its 
second ionization energy is about 18eV. Near chemical potential the conduction electron is 
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estimated to encounter a repulsive local potential Uq ~ 18. 9eVil when it scatters with the 
Zn impurity. This is much larger than the bandwidth (2eV) of the d x 2_ y 2 band of 3d Cu 2+ 
electrons. Thus, the nonmagnetic impurity Zn can be described roughly by a spin vacancy 
in the unitary limit. We start from the Hamiltonian, 



H = ~ t E PG{clc ja + h.c.)P G + Jj^(S i -S i -~n i n j )+J2Wi,i-v) 

<ij>,o <ij> i 



where I labels the site of the impurity. In the standard notation, the < ij > means the 
summation over nearest neighbors and Pq = IL(1 — n i] n ii) is the Gutzwiller's projection 
operator that prohibits double occupancy. Within the mean field approximation, the BdG 
equation is derived 
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Here u™ and v™ are the Bogoliubov amplitudes corresponding to the eigenvalue E m ; Xij and 
are the bond and resonating-valence-bond (RVB) order parameters defined by Xij — 
J2 a < c L c jo- > and Ay =< QjCjj — Ci^Cji >, respectively; 5 is the hole density. They are 
determined self-consistentl y by 
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The solution found at zero temperature had already been shown by several groupsOSE3 to 
have a nearly-zero-energy resonance for the LDOS when Uq is very large compared to J 
or t. The order parameters Ay near the impurity are suppressed and a small component 
of s-wave pairing is induced. In the slave-boson mean field theoryil the magnetic corre- 
lation obtained is overestimated. The simplest way to correct this deficiency is to use the 
eigenvectors obtained by BdG equations to construct a variational wave function with the 
projection operators rigorously imposed. For the uniform caseS a similar method has been 
used successfully. 

Following the work by Yokoyama and Shiball and Himeda et a/.@, we write this trial 
wave function for the ground state in terms of a Slater determinant of N e /2 dimension, 

l0>=^(E(^ 1 ^)^Vl i ) 7Ve/2 |O>, (8) 

ij 

where U and V in Eq.([8D are the matrices of u™ and v™, respectively. Without the Gutzwiller 
projection operator, this wave function is exactly the same as BdG ground state but with 
fixed N e electrons. The relation of this wave function with superconductivity in the absence 
of impurity was discussed in Refs. (22,24). Most properties calculated with or without the 
Gutzwiller projection operator are quite similar as shown by Zhang et a/.i3. However, the 
spin-spin correlation calculated by BdG (Eq.(2)) is very much smaller than by Eq.(f|). 

It should be noted that the trial wave function in Eq.@ is a paramagnetic RVB state 
without the antiferromagnetic long range order (AF LRO). In a uniform system without 
impurity at low doping, 5 < 0.06, this state is unstable^! with respect to the AF LRO. To take 
into account AF LRO, we could either add a Jastrow factorl!, such as exp(— h u J2i{~^-Y S l z ) , 
to modify the trial wave function or we could include spin density wave order parameteil^E 
to the original BdG equations. Since both approaches obtain almost identical results, we 
shall use a Jastrow factor here. 

In addition to the issue of AF LRO at low doping, we are also concerned with the 
lack of consideration of strong correlation in the mean field theory of BdG equations. 
Use of Gutzwiller approximation^ in BdG would improve but it still may not be enough. 



When the no-double-occupancy constraint is included exactly, we could examine the issue of 
attraction!! or repulsionEHil of holes by the impurity more accurately. Hence we introduce 
a Jastrow factor to reflect the influence of the impurity on the near-by hole distribution and 
magnetic polarization. This new trial wave function is 

>= exp (- £((-1)^ + X{1 ~ ni) ))\<t> >, (9) 
i 1 

where Ri = \J (x; — xj) 2 + (yi — yi) 2 is the distance from the impurity site denoted by /. 
The first term in the exponent in Eq.(H) introduces a spatial dependent staggered magnetic 
field, which consists of two terms, hi = h u + j*-. h u provides a uniform AF LRO at low 
doping with or without the impurity, ho is used to describe the enhanced AF correlation 
effect around the impurity. This enhancement will repel holes away from the impurity. 
Hence we include the second term associated with A for this repulsion. Notice that if A is 
negative, then the hole is attracted to the impurity and the electron is repelled from it. The 
values of h u , h and A are determined by minimizing the variational energy. In Eq.(||) we 
have chosen form to simulate the extent of the spin polarization around the impurity. We 
have examined several other functional forms and results are about the same as long as it 
covers a substantial region around the impurity. Below we will report mostly the results 
obtained from \ipj > with form. Then we will also show that similar results are obtained 
with a different trial function \ipj > using form. The latter has been previously shown 
by G. Khaliullin et a/.@ to be the spatial distribution of the impurity- induced moment. 

Our attention is also focued on the spatial magnetic polarization near the impurity. 
Without loss of generality the impurity is supposed to be situated at the center of the 
lattice. Thus, we can use the periodic boundary condition for the numerical calculation. 
For the 8x8, 12x12 and 16x16 lattice sizes we find that the spin cloud induced by the 
impurity extends only several lattice sites and all the quantities we are concerned with, 
including the local magnetization and the spin-spin corrlation function, have no qualitative 
and significant changes with the change of the lattice size. This is because the lattice sizes 
we used are large enough for the polarized spin cloud. Here we present the numerical results 
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obtained for a 12 x 12 lattice in the zero temperature limit with t/J = 3 and Uq = 100 J. 
In this paper J is our basic energy unit. We solve self-consistently the BdG equations 
and obtain the order parameters Xiji ^ij an d the BdG amplitudes U and V. The pairing 
order parameters Ay can be decomposed into extended s-wave and d-wave components as 
A d (z) = l{A x {i)+A_ x {i)-A y {i)-A_ y {i)) and A s {i) = l{A x (i) + A„ x {i) + A y {{) + A_„(z)) . 
The d-wave component is suppressed around the impurity site and it induces a small s- 
wave pairing component which is consistent with other group's results02ll. In principle, the 
Jastrow factor introduced could modify the distribution of the order parameters. In practice, 
tuning the values of the order parameters around the solutions of the self-consistent BdG 
equations has little effect on the physical quantities discussed below, except a slightly lower 
ground state energy is obtained. After obtaining the BdG solution and matrices U and 
V, we carry out the VMC simulation to determine the optimized ground state energy. 10 5 
samples were used in each MC simulation to measure the physical quantities. Since there 
are three variational parameters: h u , h and A, the calculation to find the optimal solution is 
quite involved. Here we only report the main results. In Fig.l we show energy per site as a 
function of A for two doping concentrations. In Fig.l(a) for doping concentration 8 = 0.055 
and h u = 0.05, results for ho = 0.2 (solid circles) and ho = (open cirlces) are compared. 
Fig. 1(b) shows that the lowest energy for 5 = 0.152 is achieved for h u = 0.0, ho = 0.0 and 
A = 0.2. It is noted that the energy is quite sensitive to the value of A. The lowest energy is 
acquired for positive A, thus the hole is repelled away from the impurity while the moment 
is binded to the impurity. 

We compare the optimal ground-state energy per site calculated from the trial wave 
functions \cj) > and > in Table I. In the third row we also list the total energy difference 
(AE) between these two wave functions. The variational parameters for the optimized wave 
function are listed in the table II. 

Table I: Optimal groud state energy per site as a function of hole density for two trial wave 
functions: |0 > and >. The third row lists their total energy difference 
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Doping 5 

\<P> 

IV/ > 

AE 



0.028 
-1.207±0.002 
-1.224±0.001 
-2.4±0.4 



0.055 
-1.336±0.004 
-1.351±0.001 
-2.3±0.7 



0.083 
-1.475±0.002 
-1.487±0.002 
-1.7±0.6 



0.111 
-1.611±0.003 
-1.623±0.003 
-1.7±0.9 



0.139 
-1.738±0.003 
-1.748±0.002 
-1.4±0.7 



0.152 
-1.801±0.002 
-1.806± 0.001 
-0.7± 0.3 



Table II: Optimized variational parameters for \if)j >. The values in the parenthesis are 
for the clean system without impurity. 

Doping 5 0.028 0.055 0.083 0.111 0.139 0.152 
h u 0.1(0.1) 0.05(0.05) 0(0) 0(0) 0(0) 0(0) 

h 0.3(0) 0.2(0) 0.05(0) 0.03(0) 0(0) 0(0) 

A 0.8(0) 0.6(0) 0.4(0) 0.4(0) 0.2(0) 0.2(0) 



As shown in Table I the Jastrow factor which simulates the magnetic polarization around 
the impurity in Eq.@ reduces the energy of the projected BdG wave function \<f> > by a 
significant amount. Although the energy per site has been improved only by a very small 
amount, the total energy gain is greater than 0.7J. This is a very large energy gain due 
to the influence of a single impurity. It also clearly demonstrates that BdG approach has 
significantly underestimated the magnetic correlation surrounding the impurity. 

Table II shows that h u is zero, i.e. there is no AF LRO for doping greater than 0.08 
with or without the impurity. This is expected as a single impurity cannot induce LRO for 
the whole system. At the underdoped region, for 5 = 0.083 ~ 0.11, although there is no 
uniform AF LRO, the spins around the impurity tend to form a local AF cloud as reflected 
by the nonvanishing parameter ho. It should be cautioned that in this case our trial function 
> has broken the spin up-down symmetry. A more accurate description of this state 
should be a state with a fluctuating local AF polarization but without fixing the moment 



8 



in a particular direction. For S > 0.139 although ho = and there is no apparent magnetic 
polarization around the impurity, the holes are still repelled from the impurity. This result 
disagrees with the result reported by Tsuchiura et al^. 

To examine the magnetic polarization induced around the impurity more closely, we have 
calculated the difference of the local magnetization < S Z (R) > and the spin-spin correlation 
function < S z {n)S z {n + R) > between systems with and without impurity. Both results are 
plotted in Fig. 2 as a function of the square of the distance from the impurity for several 
dopant densities. (— l) R (< S Z (R) > — < S Z (R) >o) shown in Fig. 2(a) indicates that 
< S Z (R) > is enhanced near the impurity. < .. >o is for the clean system without impurity. 
For 5 > 0.083 there is no AF LRO and the induced magnetization only exists within a few 
lattice constants around the impurity. In Fig. 2(b) we show that the spin-spin correlation 
is also enhanced near the impurity. Site n is one of the nearest neighbors of the impurity. 
Again the enhancement is weaker when the doping increases. This is consistent with the 
experimental observation. 

In Fig. (3) we plot the impurity induced spin and charge profiles, AS* 2 =< Sf > — < 
Sf > , and AN h =< > — < iV/ 1 > respectively, for two different dopant concentrations. 
Here iV/ 1 = 1 — — rij_ CT . It can be seen that the holes are kept away from the impurity 
and a spin cloud is formed around the impurity. As the hole doping increases the spin cloud 
becomes smaller in size. 

To estimate the size of the induced magnetic polarization and the induced moment, we 
calculate M(R) = 3g < \J {Y^ R {—^ySf) 2 >, where the Lande g factor g = 2 and N R is 
the number of sites within radius R of the impurity. The difference between the induced 
magnetization with and without impurity, M(R) — M (R), is plotted as a function of R 2 
in Fig.4. Results obtained by using \<f> > and > are shown in Fig. 4(a) and Fig. 4(b), 
respectively. In the inset of Fig. 4(b), results for 5 = 0.111, 0.139 and 0.152 are shown with 
a different scale. The saturation of values of M(R) — M (R) at large R indicates that the 
induced magnetization has a finite extent. We shall define the size of the induced magnetic 
polarization to be R c . At R = R c M(R) — M (R) reaches about 70% of its saturated values. 
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This moment is much larger for > than for \(f) >. Hence the local staggered magnetic 
field, hi, and the repulsion between impurity and hole introduced by the Jastrow factor in 
Eq. (|9|) has enhanced the induced moment. 

In Fig. 4(b) we have used hi = h u + ^ and for the repulsion between hole and 
impurity in Eq.(^) for >■ To examine the sensitivity of the result to the choice of the 
R dependence, we change to for both hi and the repulsion term in the Jastrow factor. 
The optimized variational energies are almost the same as the results reported in Table I. 
The results for the induced magnetic polarization is plotted in Fig. 4(c) which are quite 
similar to Fig. 4(b). 

Results in Fig. (4) show that the in the AF LRO states or 5 < 0.083, the induced magne- 
tization is much larger. When there is no LRO the induced magnetization decreases rapidly 
with increasing hole concentraion. This is consistent with experiments^. It is also consistent 
with the theoretical result reported by Tsuchiura et alM. But we do not agree with their 
conclusion that the holes are attracted toward the impurity. On the contrary, we have shown 
above that the holes are repelled away from the impurity to lower their kinetic energy. This 
effect might give an explanation to the similarity^ between Li + and Zn 2+ . The holes are also 
repelled away from the Li to gain energy. 

The induced moment M = (M(R C ) — Mq{R c ))/Nr c and the square of the size of the 
induced cloud, R\ are plotted as a function of hole concentration in Fig. 5 and its inset, 
respectively. For the hole concentration 5 < 0.055, the local magnetic moment we obtained is 
about 0.5/i£ as compared with the experimental value 0.4 ~ 1 /ig for the Zn 4% substitution 
and different dopings. The rapid decrease of the size of the induced spin cloud could be due 
to the screening by the conducting carrier 

In summary, the magnetic polarization induced by nonmagnetic impurities in high Tc 
cuprate compounds is studied by combining the variational Monte Carlo simulation and 
Bogoliubov de Gennes mean field Hamiltonian for the two-dimensional t — J model. A 
Jastrow factor is introduced to account for the induced magnetic moment and the repulsion 
between holes and the impurity. A substantial energy gain is obtained when the holes are 

10 



repelled and the antiferromagnetic polarization is enhanced near the impurity. The doping 
dependence for the induced magnetic moment is consistent with experiments. 
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FIGURES 

FIG. 1. Variational energies plotted as a function of the imurity-hole repulsion parameter A. 
Fig. (a) and (b) are for different values of the parameters as speicifed in the figure. 

FIG. 2. Enhancement of (a) S z and (b) spin-spin correlation function for diffferent hole densities 
plotted as a function of the square of the distance from the impurity. 

FIG. 3. Spin and charge profiles for different hole densities calculated by >. The parameters 
are listed in the table II. 

FIG. 4. The induced magnetic polarization, M(R) — Mq(R), plotted as a function of R 2 for 
different hole densities, obtained from (a) the BdG wave function \<p > and (b) the wave function 
\ipi >. The parameters are listed in the table II. In the inset, 5 = 0.111, 0.139 and 0.152 are shown 
with a different scale, (c) is obtained from > which is similar to > used in (b) but with a 
different functional form for the Jastrow factor as discussed in the text. The parameters used are 
5 = 0.028, h u = 0.1, ho = 0.2, A = 1.2; 5 = 0.055, h u = 0.05, h = 0.1, A = 1.0; 6 = 0.083, h u = 0, 
h = 0.025, A = 0.8. 

FIG. 5. The induced magnetic moment plotted as a function of hole density. In the inset, the 
size of the spin cloud versus hole density. The parameters are listed in Table II. 



15 



-1 .336 

-1.338 

-1 .340 

-1 .342 

P -1.344 
CD 

LLI -1 .346 

-1 .348 

-1 .350 

-1.352 - 
-1.798 

-1.799 - 

-1.800 

-1.801 

D) -1.802 



CD 
LU 



-1 .803 
-1.804 - 
-1.805 
-1.806 
-1 .807 



h u=0.05, ho=0.2 
h u =0.05, h =0 

5=0.055 



(a) 




h,,_ 



u=0, ho=0, 5=0.152 



(b) 




-0.2 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



1.2 



16 




17 




18 



25 

20 

15 

10 

5 



-5 

10 
70 

60 

50 

40 

30 

20 

10 


60 

50 

40 

30 

20 

10 



10 





20 



